Unconventional anisotropic superexchange in a'-NaV 2 05 
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The strong line broadening observed in electron spin resonance on NaV20s is found to originate 
from an unusual type of the symmetric anisotropic exchange interaction with simultaneous spin- 
orbit coupling on both sites. The microscopically derived anisotropic exchange constant is almost 
two orders of magnitude larger than the one obtained from conventional estimations. Based on this 
result we systematically evaluate the anisotropy of the ESR linewidth in terms of the symmetric 
anisotropic exchange, only, and we find microscopic evidence for precursor effects of the charge 
ordering already below 150 K. 
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The isotropic exchange constants in one-dimensional 
antifcrromagnets are obtained by measurements of the 
magnetic susceptibility or by inelastic neutron scatter- 
ing. The anisotropic exchange contributions, however, 
are only accessible by means of electron spin resonance 
(ESR) , because the spin-spin relaxation measured by the 
ESR linewidth is driven primarily by the correspond- 
ing effective local fields. Conventional theoretical esti- 
mations of anisotropic exchange parameters yield values 
by far too small to describe the experimental results. A 
prominent example for this problem is the spin-ladder 
a'-NaV^Cv This compound was initially identified as a 
spin-Peierls system 0] which triggered enormous efforts 
to investigate the nature of this transition, ft was found 
that the system actually undergoes a charge-order (CO) 
transition at Tqo ~ 34 K Q from a uniform oxidation 
state of V A - b+ ions at high temperature Jj into a state 
with "zig-zag" type charge distribution |j] accompanied 
by spin-gap formation. Moreover, various experimental 
studies revealed an anomalous behavior at about 200 K, 
far above Tco , which has been attributed to the existence 
of charge fluctuations in the system 0, 0, • 

Previously, it was proposed that the spin relaxation 
in NaV2C>5 is strongly affected by these charge fluctua- 
tions 8]. ESR directly probes the spin of interest and, 
hence, is extremely sensitive to such dynamic processes 
of the electronic structure. The underlying mechanism 
of the spin relaxation, however, still remained a matter 
of heavy debate |J llflj . fn this work, we identify the 
anisotropic exchange (AE) interaction as the dominant 
source of line broadening and we calculate the AE pa- 
rameters on the basis of a microscopic charge-distribution 
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picture. With these parameters we are able to describe 
the angular dependence of the ESR linewidth AH at tem- 
peratures above Tco- The resulting temperature depen- 
dence of the exchange parameters is a clear fingerprint of 
the increasing charge fluctuations on approaching Tco- 

All details concerning the preparation and characteri- 
zation of the crystals and the experimental ESR set-up 
have been published previously [llj. The observed ESR 
signal in a-NaV^Os consists of a Lorentzian line with a 
g-value g rs 2, characteristic of a spin-only system with 
quenched orbital moments ^3- The linewidth increases 
monotonously from a value of 10 Oe at Tco up to sev- 
eral hundredth Oe above room temperature, with the 
linewidth AH C for the magnetic field applied along the 
crystallographic c axis being about twice as large as AH a 
and AHjj (along the a and b axis). 

At first, let us consider possible origins for the line 
broadening AH in NaV20s. Single-ion (5=1/2), hyper- 
fine and spin-lattice relaxation were shown to be less 
important in NaV 2 5 (0 El E3) The anisotropic 
Zeeman-effect is not relevant, because of nearly equiv- 
alent g tensors for all vanadium sites. Therefore, only 
three sources remain to account for the broadening of 
the ESR spectra in NaV^Os - the dipole-dipole (DD), the 
symmetric anisotropic-exchange (AE) and the antisym- 
metric Dzyaloshinsky-Moriya exchange (DM) interac- 
tions. These contributions have been already estimated 
and discussed in Ref. The assumption that the DM 
interaction is the main perturbation was based on the 
conventional relation for the DM-vector \d\ ~ {Ag/g)\J\ 
(where g and Ag are the g factor and its anisotropy, re- 
spectively) ^3 • The ESR data could be modeled by this 
approach, but it was necessary to assume the presence 
of strong charge disproportions even at highest temper- 
atures to allow for an appropriate direction of the DM 
vector 8]. Later, Choukroun et at Q questioned the 
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FIG. 1: Possible paths for AE between two sites A (with 
ground state jjy), exited state £)) and A (|£) and \cp), re- 
spectively). Solid arrows correspond to the effective hopping 
integrals, dashed arrows indicate the matrix elements of the 
spin-orbit coupling. Numeration corresponds to the sequence 
of the matrix elements in forth order perturbation expansion 
(see Eq.QJ. 




dominance of the DM interaction, showing that the con- 
tribution of the DM interaction to the ESR linewidth 
in quantum-spin chains cannot be larger than that of the 
AE. But the AE itself as taken from conventional estima- 
tions is by far too small to account for the large linewidth 
observed in NaV2C>5. Therefore, such conventional es- 
timations have to be taken with care. A recent field- 
theoretical treatment of quasi-one-dimensional S = 
antiferromagnetic chains came to similar conclusions 
because the DM interaction was found to produce a di- 
vergence in the temperature dependence of the linewidth 
AFdm ~ T~ 2 for T < J/k B - This is in contrast to the 
monotonic increase of AH with increasing temperature 
in NaV2C>5 |2J. Such a behavior, however, is in agreement 
with the theoretical expectation for a dominant AE . 
Experimental investigations of related compounds cor- 
roborate this expectation, too 0, HH Il7l lis| . In this 
respect LiCuV04 received a key role because the 
DM interaction can be completely ruled out by its crys- 
tal symmetry. The linewidth is dominated by the AE, 
because the ring-exchange geometry in the CU-O2 chains 
strongly enhances the AE as compared to the conven- 
tional estimation |Jae| ~ (Ag/g) 2 \J\ Q. In the fol- 
lowing we will provide detailed microscopical estimations 
of this term in NaV2C>5 and show that the angular and 
temperature dependencies of AH can be completely de- 
scribed in terms of this relaxation mechanism, only. 

Starting with the microscopic analysis of the AE paths, 
the Hamiltonian for AE between two neighbouring sites 
A and A' can be written as Hae — S% ■ D^p • S^, , where 
{ev,/3} = {x,y,z}. Taking into account all possible vir- 
tual processes (displayed schematically in Fig.[TJ between 



FIG. 2: Schematic pathway of AE between V ions in NaV^Os. 
Big spheres denote V ions, small spheres O ions. Frame (a): 
intra-ladder exchange between the ground state d xy orbitals. 
Frame (b): intra-ladder exchange between the excited d x 2_ y 2 
states. Frame (c): Possible inter-ladder exchange paths. 



site A to site A', we derive the following expression for 
D^g in forth order of perturbation theory: 
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where, e.g., is the effective hopping integral between 
the states |£) and |C) via intermediate oxygens, and 
(£,\l a \() denotes the matrix element of the spin-orbit (SO) 
coupling Hso = A/ Q s a . Here we assume that the charge- 
transfer energy Aaa' from site A to site A' is large com- 
pared to the crystal-field splittings A c f = A^ V ,A^ V . 
The first two correspond to conventional AE processes 
[H El, while the others (Fig. Hfc-f)) and the general 
expression Q are presented to the best of our knowledge 
for the first time. For example, in case (f ) the electron at 
site A is excited via SO coupling from the ground state 
I rf) into the state \C), then it is transferred to the empty 
state \(p) at site A' and interacts via SO coupling with the 
electron in the corresponding ground state Finally, 
one of the electrons hops from state |£) to the initial state 

\v)- 

Focusing now on NaV205, we recall that the elec- 
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FIG. 3: Angular dependence of the ESR-linewidth at different 
temperatures. Fit curves (lines) are described in the text. 
Upper frame: normalized to the linewidth for the magnetic 
field applied along the c axis. Lower frame: illustration of the 
contributions of intra- and inter-ladder AE to the linewidth 
far above (dashed line) and near Tco (solid line). 



tron is distributed between two V-ions on the same 
rung. Correspondingly, its ground state wave func- 
tion 1 77) is given as a superposition c\ \d xy ) — c 2 \d xy ) of 
the two vanadium (i-orbitals. Analogously, the ground 
state of the electron on the adjacent rung is given by 
|£) = c[ \d xy ) — c' 2 \d xy ). We illustrated the correspond- 
ing ci-orbitals together with relevant bridging oxygen p- 
orbitals (7r-bonding with hopping integral t^ v = t n ) in 
Fig. Ufa) for the high-temperature limit, where all coef- 
ficients ci,C2,c[,C2 become equal to l/v2- Note, that 
due to the orthogonality of the wave functions processes 
(a)-(d) (Fig. [TJ do not contribute to AE within one lad- 
der in NaV2C>5. Therefore, we will now concentrate on 
processes (e) and (f) and discuss the relevant excited 
states |C) and \<p) involved. Considering the possible ex- 
citations of the electrons via SO coupling we find that 
the largest contribution is obtained by the matrix ele- 
ment (d xy \l z \d x 2_ y 2) — 2i. Hence, the relevant excited 
states are the combinations c\ \d x 2_ y 2) — c 2 \d x 2_ y 2) and 
d x \d x 2_ y 2) — c' 2 \d x 2_ y 2) for A and A' rungs, respectively. 
The charge-distribution picture for the excited states 
(cr-bonding via oxygen p-orbitals with hopping integral 



t' a ) is shown in Fig.^b). 



Thus, using expression one can derive D zz as 
t t' 
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To estimate D zz we use the free ion value A = 31 meV 
|2l|. the splitting between the d xy and the d x 2_ y 2 states 
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FIG. 4: Right ordinate: temperature dependence of the 
linewidth-ratio for the magnetic field applied along a or b axis 
normalized to AH C . Left ordinate: temperature dependence 
ratio of the inter- to intra-ladder AE constants obtained from 
the fitting of the angular dependencies of AH . 
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0.36 eV l22l | 2 3l. and the charge-transfer en- 
3 eV [24(. The hopping integral t' a can- 
not easily be calculated, however, one can assume t' a w 
<7r = 0.17 eV Q as a lower bound for t' a , and we ob- 
tain D zz « 0.6 meV in the high-temperature limit where 
the electron is equally distributed on each rung. This 
yields a characteristic linewidth AH ~ 300 Oe in very 
good agreement with the experimental linewidth. Note 
that our estimate is about two orders of magnitude larger 
than previous results |l2l | , corroborating the importance 
of microscopic considerations for the estimation of AE 
parameters. 

Taking now into account possible inter-ladder ex- 
change paths as shown in Fig. Efc), we obtain contri- 
butions of comparable strength for the inter-ladder AE. 
These paths involve a 90°-exchange geometry, which has 
been discussed in detail in Refs. Hil I25I Etl Considering 
all possible inter-ladder exch ange paths in the appropri- 
ate local coordinate systems [17J, |2jj we find a maximal 
component of the effective AE tensor along the crystal- 
lographic a axis. 

Having identified and estimated the source of the ESR 
line broadening in NaV20s, we will now apply this model 
to the experimental data. The analysis of the angu- 
lar dependencies in terms of second moments has been 
discribed previously |l7t Il8j . The experimental angular 
dependencies [5^ of AH are shown in Fig. |3 together 
with the fit curves. As a result we derive the ratio of 
the two essential fit parameters -Di n t ra and -Dintm- lor 
the AE parameters within and between the ladders, re- 
spectively. Fig. 0] shows the temperature dependence of 
the ratio Di nter / -Dintra together with the linewidth ratios 
AH a /AH c , AH b /AH c (note that only the ratio of the 
exchange parameters can be determined from the ESR 
linewidth at temperatures T < J/ks as discussed in 
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Ref. Il7l) . It can be clearly seen that at high temper- 
atures (T > 150 K) the dominant contribution to the 
line broadening is given by the intra-ladder AE. On de- 
creasing temperature, below 150 K, the ratio strongly 
increases and the inter-ladder contribution becomes dom- 
inant. This can be understood taking into account the 
strong dependence of the AE parameters from the co- 
efficients ci , d x , C2 , c' 2 which describe the electronic oc- 
cupation on the vanadium sites. That means e.g. for 



D zz (Eq. EJ) the coefficient 



2 L 2. 



is equal 



to 1 for the case F 4 - 5+ — V 4 ' 5+ , and vanishes for the "zig- 
zag" charge order (V 5+ — V 4+ ) realized below Tqo |4J. 
The observed increase of the -Dintcr/Mntra ratio already 
far above Tco indicates that precursors of the develop- 
ing CO set in at about 150 K, weakening considerably 
the intra-ladder AE. 

Further evidence for the onset of charge disproportions 
above Tco has been reported by the strong frequency de- 
pendence of the ESR linewidth between 34-100 K |2?Sj . 
and the anomalous features observed by optical spec- 
troscopy measurements 0, 0, H3| • Note that the uniform 
susceptibility has significant deviations from the Bonner- 
Fisher law already at T < 200 K, too |7j. The coupling of 
these charge fluctuations to the lattice has been revealed 
by a softening of the elastic constants below 100 K de- 
tected by ultra-sound experiments 0, and by the 
shift of the phonon energy found in light-scattering mea- 
surements around 80 K pBSf. Moreover, we would like 
to point out that similar observations in CuGeOs have 
been explained in terms of lattice fluctuations existing 
already far above the spin-Pcicrls transition |l7j . We be- 
lieve that the proposed spin-relaxation mechanism and 
the microscopic estimations do not only apply for the 
case of NaV2 05, but may allow to describe the spin dy- 
namics in many transition-metal compounds. 

In summary, we have identified the symmetric 
anisotropic super-exchange to be the source of the im- 
mense ESR line broadening in NaA^Os. In this micro- 
scopic picture the dominant process consists of the simul- 
taneous virtual hopping of electrons between the ground 
states and excited states of vanadium ions on neighbor- 
ing rungs of the ladder involving the spin-orbit coupling 
on both rungs. This novel unconventional exchange pro- 
cess has not been considered in the discussion of ESR 
line broadening before. The corresponding AE param- 
eter is found to be of the order of 1% of the isotropic 
exchange constant resulting in a high temperature limit 
of the ESR linewidth of approximately 10 2 Oe. On the 
basis of this microscopic analysis we have shown that 
the ESR data can be entirely described in terms of the 
symmetric anisotropic exchange only. The temperature 
dependence of the linewidth and derived exchange pa- 
rameters evidences the presence of charge fluctuations in 



NaV205 up to 150 K on a microscopic level. 
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